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Abstract 

We consider a two-dimensional Coulomb gas of positive and negative pointlike 
unit charges interacting via a logarithmic potential. The density (rather than the 
charge) correlation functions are studied. In the bulk, the form-factor theory of an 
equivalent sine-Gordon model is used to determine the density correlation length. 
At the surface of a rectilinear plain wall, the universality of the asymptotic behavior 
of the density correlations is suggested. A scaling analysis implies a local form of 
the compressibility sum rule near a hard wall. A symmetry of the Coulomb system 
with respect to the Mobius conformal transformation, which induces a gravitational 
source acting on the particle density, is established. Among the consequences, a 
universal term of the finite-size expansion of the grand potential is derived exactly 
for a disk geometry of the confining domain. 
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1 Introduction and summary 

In this paper, we consider classical Coulomb systems in thermodynamic equilibrium. For 
the sake of simplicity, we will restrict ourselves to the case of a symmetric two-component 
plasma (TCP), i.e., a neutral system of two species of particles of opposite unit charges 
a — ±1, living in a //-dimensional space and interacting through the Coulomb interaction 
Oi<jjv(\vi — r j I ) . In dimension v, the Coulomb potential v at a spatial position r e R u , 
induced by a unit charge at the origin, is the solution of the Poisson equation 

Av(r) = -s v 6(t) (1.1) 

where s u is the surface area of the //-dimensional unit sphere; s 2 = 2n and s 3 = 4n. In 
particular, 

v(r) = -ln(|r|/r ), u = 2 (1.2a) 

v(r) = l/|r|, i/ = 3 (1.2b) 

where, for v — 2, the length scale r will be set to unity without any loss of general- 
ity. For the case of pointlike particles, the singularity of v(r) at the origin prevents the 
thermodynamic stability against the collapse of positive-negative pairs of charges: in two 
dimensions for small enough temperatures, in three dimensions for any temperature. 

Let us now introduce some notations. We will work in the grand canonical ensemble, 
characterized by the fixed domain D in which the TCP is confined, by the inverse tem- 
perature P, and by the couple of equal (constant) particle fugacities z + = Z- = z. The 
microscopic densities of charge and of the total particle number are defined respectively 
by 

P( r ) = Z>Mr), ™( r ) = I>-( r ) ( L3 ) 
where n a (r) = Y<i&a,<Ji&{ r — r %) is the microscopic density of particles of species a — ± 
and % indexes the particles. The thermal average will be denoted by (• • •). At one particle 
level, the total charge and particle number densities are given respectively by 

p(r) = (p(r)>, n(r) = (n(r)> (1.4) 
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Due to the charge ± symmetry, n+(r) = n_(r) = n(r)/2. At two-particle level, one 
introduces the two-body densities 

7w(r,r') = I Y^Sa^a'^K* ~ r i)$i r ' ~ r j) ) 

= (n CT (r)rv(r')) - (h a (v))5 aa ,5(r - r') (1.5) 

Clearly, n ++ (r, r') = n (r, r') and n + _(r,r') = n_ + (r, r'). The corresponding Ursell 

functions read 

U aa >(r,r') = rw(r,r') - n CT (r)n CT /(r') (1.6) 
They will occur in the charge and density combinations 

U p (r,r') = XWtf™'(r,r') (i. 7a ) 
C/ n (r,r') = J2U a(T ,(r,r') (1.7b) 

respectively. The truncated charge-charge and density- density structure functions are 
defined by 

S p (r,r') = (p(r)p(r'))-(p(r))(p(r')) (1.8a) 

S n (r,r') = (h(r)h(r'))-(h{r))(h(r')) (1.8b) 

respectively. It is useful to consider also the pair correlation functions 

^(^^/(r') 

in their charge and density combinations, defined by 

M r > r = ^E^Tw^r') ( L10a ) 



<7,<7' 



h n {v,v') = \Y, h °A*J) ( L10b ) 

As any Coulomb system, the TCP admits a description in the Debye-Hiickel (high- 
temperature) limit [I]]. 

In two dimensions and for the case of pointlike charged particles, the particle density 
n is an "irrelevant" variable which only scales the distance. The complete description of 



the bulk thermodynamics is available in the whole range of inverse temperatures (3 < 2 
where the plasma is stable against the collapse of positive-negative pairs of charges. The 
exact equation of state 

(3p = n(l- (3/4) (1.11) 

where p is the pressure and n the total particle density, has been known for a very long 
time 0, |J. The evaluation of other thermodynamic quantities (free energy, internal 
energy, specific heat, etc.) can be based on an explicit density- fugacity (n — z) rela- 
tionship. This relationship was obtained only recently |4]] via a mapping onto a classical 
two-dimensional sine-Gordon theory. Surface thermodynamics (surface tension) of the 
two-dimensional TCP in contact with a rectilinear dielectric wall was derived for specific 
boundary conditions in refs. [||, |6j. The large-distance behavior of the bulk charge-charge 
correlation function was given in ref. pj by exploring the form-factor method for the 
equivalent sine-Gordon theory. Exact formulae for pair correlations for an arbitrary in- 
terparticle distance are available just at the collapse point (3 = 2 || [J, which corresponds 
to the free-fermion point of an equivalent Thirring model. 

The long-range tail of the Coulomb force causes screening, and thus give rise to exact 
constraints for the charge correlations for any value of (3 (for a review, see ref. fToH). 

In the bulk, the charge structure function S p (r, r') = S p (\r — r'|) obeys the Stillinger- 
Lovett sum rules [JTTJ] , namely the zeroth-moment condition 



dr S p (r) = (1.12) 



and the second-moment condition 



dr \r\ 2 S p (r) = -— -, u = 2,3 (1.13) 

(3tt(p — 1) 

Hereinafter, without any loss of generality, we put the dielectric constant e of a medium 
in which the plasma lives equal to unity. 

Let us now introduce a semi-infinite TCP which occupies the half space x > 0, and 
denote by y the set of (y — 1) coordinates normal to x. The plane at x = is a hard wall 
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impenetrable to the particles. The half-space x < is assumed to be filled with a material 
of dielectric constant ew'- a particle of unit charge at the point r = (x > 0, y) has an 
electric image of charge (1 — ew)/(l+ew) a ^ the point r* = {—x, y) £L2|. Due to invariance 
with respect to translations along the wall and rotations around the x direction, 



S p (r, r') = S p (x, x'\ |y - y'|) = S p (x', x; |y - y'|) 



;i.i4) 



The electroneutrality condition Ql .12 ) takes the form 

/"CO r 

J dx' J dy S p {x,x';y) = (1.15) 
The Carnie and Chan generalization of the second-moment condition ( TTT^ ) results 



in the dipole sum rule M, [15 



J dx J dx' J dy (x' — x)S p (x, x'; y) 



v = 2,3 



;i.i6) 



2Pn(u - 1)' 

The asymmetry of the screening cloud of a charged particle sitting near the wall induces a 
long-range tail in the charge correlation along the wall [[U| [17) , except of the cases ew = 
(ideal dielectric wall) and ew = oo (ideal conductor wall) which will be excluded from our 
considerations. One expects an asymptotic power-law decay 



S p (x,x';y) 



/p(x, x ) 



CX) 



;i.i7) 



where f p (x,x') obeys the sum rule |L8], [19 



dx / dx f p (x, x ) — , . .vin) 
is jo 2p[7r(i/ - 



u = 2,3 



;i.i8) 



Recently [Pm , it was proven that for any value of x > it holds 
dx' / dy (x'-x)S , p (x,x / ;y) = — ir{v - 1) / dx'f p (x,x'), 



v = 2,3 



;i.i9) 



When both sides of ( 1.19j) are integrated over x from to oo, it is clear that the sum rule 
( |1.18| ) for f p is a direct consequence of the dipole sum rule ( |1.16| ), and vice versa. 

As concerns the density correlation function, according to the general theory of fluids 



21| , the zeroth moment of the bulk density structure function S n (r, r') = 5 n (|r — r'|) is 



related to the isothermal compressibility 

dp) 



1 (dn\ 



n 



via 



h I dr Sn{r) = \ xp (L2i: 



In two dimensions and for the case of pointlike charged particles, with the use of the 
exact equation of state ( |1.11| ), one gets explicitly |22| 



d2r s - w = T^rm (1 - 22) 



Recently, using the technique of a renormalized Mayer expansion [p3|| , and in particular a 
remarkable "cancelation property" of specific families of renormalized diagrams [24], the 
second moment of S n was shown to have a simple value |25| 



/ dV|r|2S " ( ^lM^W (L23) 

The universal finite-size properties of two-dimensional critical systems with short-range 
interactions among constituents are well understood within the principle of conformal 



invariance []26] - [50]. For a finite system of characteristic size R, at a critical point, 
the dimensionless grand potential (3Q = — In S (S is the grand partition function) has a 
large- R expansion of the form 



/3Sl = AR 2 + BR-—\nR + const + • • • (1.24) 
6 

The coefficients A and B of the bulk and surface parts are non-universal. The coefficient 
of the lni?-term is universal, dependent only on the conformal anomaly number c of 
the critical theory and on the Euler number \ of the manifold on which the system is 
confined. In general, % = 2 — 2h — b, where h is the number of handles and b the number 
of boundaries of the manifold (x = 2 for a sphere, x = 1 f° r a disk, x — for an anulus 
or a torus). The grand potential of the two-dimensional TCP is supposed to exhibit a 
universal finite-size correction of type ( |1.24| ) at any temperature of the conducting regime. 
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Plausible, but not always rigorously justified, arguments for a critical-like behavior were 
first given for Coulomb gases with periodic boundary conditions then for Coulomb 
systems confined to a domain by plain hard walls p2| , by ideal-conductor walls [[3^] and 
finally by ideal-dielectric boundaries [35| . The explicit checks were done at the exactly 
solvable (3 = 2 inverse temperature for various geometries of confining domains. Only very 
recently , a direct derivation of the universal finite-size correction term was done for 
the specific case of the TCP living on the surface of a sphere of radius R. By combining 
the method of stereographic projection of the sphere onto an infinite plane with the linear 
response theory, the prefactor to the universal lni? correction term was related to the 
bulk second moment of the density structure function S n , eq. ( |1.23| ). The obtained 
result confirms the prediction ( |1.24|) for a Coulomb system, as if we had c = — 1, in full 



agreement with heuristic approaches and exact results at (3 = 2 [31|-[35 



In general, due to screening phenomena, the sum rules for the charge structure function 
S p are not modified by a short- distance regularization of the Coulomb potential. On the 
other hand, the moments of the density structure function S n depend on the particular 
form of the short-range particle interactions. For the case of two dimensions and the 
pointlike character of charged particles on which we shall concentrate in this work, the 
zeroth-moment ( 1.12 ) and the second-moment ( |1.13| ) Stillinger-Lovett conditions for S p 
have their exact counterparts Ql.22 ) and ( 1.23[) , respectively, for S n . Instead of screening, 
the scaling properties and the critical-like state are relevant. The aim of this paper is to 
document how these features of the two-dimensional Coulomb fluids manifest themselves 
in the density correlation functions (in the bulk or close to a boundary) and in the universal 
finite-size correction term of the grand potential. 

Section 2 is devoted to the density correlations in the bulk. The Debye-Hiickel (3 — > 
limit is presented in subsection 2.1, the exactly solvable (3 = 2 case is treated in subsection 
2.2. A general analysis, based on the form-factor theory of the equivalent sine-Gordon 
model in analogy with ref. Jj]], is given in subsection 2.3. For comparison, the correspond- 
ing formulae for the charge correlations are presented, too. 
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Density correlations near a rectilinear hard wall are analysed in section 3. In subsection 
3.1., we present the results of the Debye-Hiickel j3 — > limit. Subsection 3.2. deals with 
the exactly solvable (3 = 2 case for the plain hard wall with ew = 1. In comparison 
with formula (|1.17|) for S p taken at dimension v = 2, a more rapid, but still power-law, 
asymptotic decay of S n along the boundary is observed, 

S n (x,x';y)~l^2 (1.25) 
T 

In subsection 3.3, at least for the considered plain hard wall, the universal form of the 
function f n (x, x') at the boundary x = x' = is conjectured, 

fn(0,0) = 7^- 2 (1-26) 

Here, we also suggest explicit forms of f p (x,x') and f n (x,x') for any value of (5. 

In section 4, we explore the scaling properties of the two-dimensional TCP confined to 
a disk of radius R, for the sake of simplicity by an uncharged plain hard wall. We present 
some important formulae which are used in the subsequent sections. The rectilinear wall 
is obtained as the limiting R — > oo case of the disk. As a by-product of the formalism, 
we derive a local form of the compressibility (zeroth-moment) sum rule ( 1.22 ) for S n near 
the rectilinear hard wall. 

In section 5, by using a specific (namely Mobius) conformal transformation of particle 
coordinates we show how the two-dimensional TCP can be mapped onto the one under the 
action of a gravitational source, acting in the same way on both positively and negatively 
charged particles. The mapping implies a new sum rule for S n in the disk geometry. In 
this section, the new sum rule is used to derive the counterpart of the bulk second-moment 
sum rule for S n ( 1.23| ) near a rectilinear wall. 



Section 6 deals with the universal finite-size correction In i?-term for the disk geometry. 
Using the new sum rule for S n derived in section 5, the universal term is confirmed at any 



value of p. This term has already been calculated at P = 2 in ref. [32 



A brief recapitulation and some concluding remarks are given in section 7. 



8 



2 Bulk density correlations 



2.1 Weak coupling 

In this subsection, we derive the asymptotic form of the density correlation function h n (r) 
in the bulk, at the lowest order in (3. It should be noted that /3-expansions of correlation 
functions must be taken for a fixed value of the inverse Debye length k = (27r(3n) 1 ^ 2 , since 
k only fixes the length scale. 

In the present case of a charge-symmetrical two-component plasma, the ordinary 
Ornstein-Zernike (OZ) equation split into two independent relations for the charge-charge 
and density-density functions |25| 



hp = Cp + c p *n*h p (2.1a) 
h n = c n + c n *n*h n (2.1b) 

respectively, where * denotes a convolution product. In the renormalized Mayer expansion 
@, |25|] , at the lowest order in /3, the charge direct correlation function is given by 

c,(l,2)= 0----0 =-/Ml,2) (2.2) 

with v being the Coulomb potential, and the density direct correlation function is given 
by the Meeron diagram 

c ^ 2 )= fZ> 2 =k K2 ^ 2) (2 - 3) 

The wavy line denotes the renormalized bond K (sum of chains). In the bulk, K(r) = 
—/3K [kt), where K is a modified Bessel function. 

Inserting (|2.2j ) into (|2.1a ), one gets the definition of K. This is why 



h p (r) = -(3K (nr) ~ -(3 ^ exp(-Kr) (2.4) 

since Kq{x) has the asymptotic form [n / '(2a;)] 1 / 2 exp(— x). In the OZ relation (|2.1b|) , where 
c n ( p.3|) is of order f3 2 , the convolution term is easily seen to be of higher order j3 3 . Thus, 
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at lowest order in (3, h n {r) = c n (r), 



j3 2 nf3 2 
M r ) = —KI(kt) ~ -— exp(-2«r) 

The same analysis can be done in dimension v = 3, where 



K{r) = — — exp(— pzr) 
r 



(2.5) 



(2.6) 



where now K = (4 7 r/3n) 1 / 2 [[37]]. In this case, the large-r behavior is not touched by an 
inevitable short-distance regularization of the Coulomb potential. 



2.2 (3=2 

When {3 = 2, for a fixed fugacity z, the particle density n — > oo and {h p , h n } — > 0. 
However, the Ursell functions U p = n 2 h p and U n = n 2 h n are exactly known [^, [J as 



U P (r) 



and 



m 
2^ 



m 



2\ 2 



K 2 (mr) + K^imr) 



U n (r) =2(^-1 [^(mr) - K 2 (mr) 



(2.7) 



(2i 



where m = 2irz (z is the fugacity). Replacing the modified Bessel functions by their 
asymptotic expansions 



K (x) = 
K x {x) = 

gives the asymtotic form of U p (r) as 



U P (r) 



and the asymtotic form of U n {r) as 

U n (r) 



2x/ 





1 




1 - 




+ ... 




8x 




1 + 


3 






+ ... 




8x 





(2.9a) 
(2.9b) 



m 



2nr 



exp(— 2mr) 



(2.10) 



m 

4 7rr 2 



exp(— 2mr) 



(2.11) 
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2.3 Form-factor analysis 

We now determine the large- distance behavior of h n (r) in the whole stability interval of 
inverse temperatures < (3 < 2 via a form-factor analysis. 



The grand partition function of the TCP can be turned into (see, e.g., [38 

^_ exp (-£(*)) 

" fV<f> exp (-S'(O)) 

where 



(2.12) 



= / d 2 r 



' (V0) 2 -2zcos(60) 



(2.13a) 



.167T 

6 2 = (3/4 (2.13b) 

is the Euclidean action of the classical sine-Gordon model. In the sine-Gordon represen- 
tation, the density of particles of one sign a = ± is 

n a = z a {e 1 ^) (2.14) 



where (■ ■ •) denotes the averaging over the sine-Gordon action (|2.13 ), two-body densities 
7SD are expressible as follows 

V( r .O = z^ie^e^') (2.15) 

etc. The parameter z = z + = z_, which is the fugacity renormalized by a (diverging) 
self-energy term, gets a precise meaning when one fixes the normalization of the cos-field. 
For the TCP Q, this normalization is given by the short-distance behavior 

n + _(r, r') ~ z + z_\r — r'|~^ as |r — r'| — > (2-16) 

dominated by the Boltzmann factor of the Coulomb potential. The corresponding formula 
in the sine- Gordon picture 

(e^«e- i6 ^ r ') ~ |r - r'l^ 2 as |r - r'| -»■ (2.17) 

is known in quantum field theory as the conformal normalization. 
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The sine-Gordon model (|2.13|) is massive in the region < b 2 < 1 (0 < < 4). It 



is integrable |3J|: its particle spectrum consists of one soliton-antisoliton pair of equal 
masses M and of soliton-antisoliton bound states ("breathers") {Bj;j = 1,2, ... < l/£}. 
Their number depends on the inverse of the parameter 

b 2 ( (3 \ 



q 1-b 2 
The mass of the .Bj-breather is given by 



;2.i8) 



rrij = 2Msin — j 



;2.i9) 



and the breather disappears from the spectrum just when rrij = 2M. Under the conformal 
normalization (|2.17|) , the relationship between the soliton mass M and the parameter z 



was established in 40 



T(6 2 



7rr(l - b 2 ) 



M 



v^r((i + o/2) 



2-2b 2 



[2.20) 



2IW2) 

where T stands for the Gamma function. Using the Thermodynamic Bethe ansatz, the 
specific quantity 



lim — In [ 



8sin(7rO 



was found in ref. ET[. As a thermodynamic result 



M 2 (nC 

n = — r^rtan — 

4(1 ~b 2 ) \2 / 



(2.21) 



[2.22) 



Note that as 3 approaches the collapse value 2, for a fixed z, M is finite and n — > oo as 
it should be. 

For the underlying sine-Gordon theory, the two-point truncated correlation functions 
of local operators O a (a is a free parameter) can be formally written as an infinite con- 
vergent series over multi-particle intermediate states, 



(O a (r)O a ,(r')) T 



AT' 

JV=1 JV • ei,...,e N 



°° de 1 ...de 



N 



(2tt 



On, 



e\...e N 



A' 



X 



£jv...ei 



a>\VN, 



0!) exp [ -|r - r'| ^ m £j cosh^- ] (2.23) 

3=1 
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where e indexes the particles [e = +(— ) for a soliton (antisoliton) and e = j for a breather 
Bj (j = 1,2,.. .)] and 9 is the particle rapidity. The form factors 

F a {9 1) ... ) 9 N ) ei ... eN = (0\O a (0)\Z ei (di),...,Z eN (9 N )) (2.24a) 
eN - 61 F al (6 N ,...,e 1 ) = (Z eif (9 N ),...,Z ei (9 1 )\O af (p)\0) (2.24b) 

are the matrix elements of the operator at the origin, between an iV-particle in-state (as 
a linear superposition of free one-particle states \Z t (6))) and the vacuum. 

In the limit |r — r'| — ► oo, the dominant contribution to the truncated correlation 
function in ( |2.23| ) comes from a multi-particle state with the minimum value of the total 
particle mass Y,jLi m €j, at the point of vanishing rapidities 9j — > 0. Due to topological 
reasons, solitons and antisolitons coexist in pairs, the total mass of the pair being 2M. 
The breathers Bj with masses given by relation ( |2.19| ) are lighter and therefore, when 



they exist and their form-factor contributions do not vanish, they are the best candidates 
for governing the asymptotic behavior of the two-point correlation function. 

The form factors of an exponential operator O a = exp(ia0(r)) for various combinations 
of particles were calculated in refs. [42| - [p4 ]. The one-breather form factors, which do 



not depend on the rapidity, posses the following general structure 

(Ole^l^)) = (5,(0)1^10) oc a j (e ir7b *) sin {^j (2.25) 

where the full dependence on a = ±1 is presented. These form factors are nonzero only 
if j = odd integer. Due to the invariance of the sine-Gordon action (|2.13 ) with respect to 



the transformation, it holds (e lbr ^) = (e _lb< ^). With regard to (|2.15|) , the total contribution 
of a given breather Bj (j odd) to the charge h p and density h n correlation functions Q1.10Q 
is proportional to 



h p oc Yl ™'( \ ei ° H \ B M)( B MK' H \ ) ( 2 - 26a ) 

<7,er'=±l 

K oc Y (OK^lB^iB^le^lO) (2.26b) 

cr,cr'=±l 

Inserting ( |2.25| ) into ( |2.26| ) one observes that one-breather states contribute only to 
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hp. Using Q2.22j ) in (|2.19|) , the mass of the lightest (elementary) 5i-breather is 



mi = k 



sin(7r/V(4-/3)) 



1/2 



[2.27) 



tt/3/(4-/3) 

where k denotes as usually the inverse Debye length. At asymptotically large r, (|2.23|) 
gives 

h p {r) oc exp(-mir), < (3 < 2 (2.28) 



in agreement with the (3 — > limit Q2.4|) . At the free-fermion point (3 = 2, the i^x-breather 
disappears, and the soliton-antisoliton pair with mass 2M determines the correlation 
length: 

U p (r) oc exp(-2Mr), (3 = 2 (2.29) 

From ( ggg ), M = 2nz = m at (3 = 2 (b 2 = 1/2), and the asymptotic form (|2~T0|) is 
reproduced. Note that from ( |2.19|) mi — > 2M as /3 — > 2, and so the inverse correlation 
length varies continuously near (3 = 2. The explicit inverse-power law dependence of 
prefactors of asymptotic formulae ( p.28|) and (|2.29|) was presented in ref. 0. 

At small (3, the large-distance behavior of h n is determined by the two-Si-breather 
state. Indeed, the corresponding form factor 0, |43| 



\e i - b f\B 1 (e 2 ),B 1 (e 1 )) = (B 1 (e 2 ),B 1 (e 1 )\e M *\0) «a 2 (e^) 



;2.30) 



(where the full dependence on a = ±1 is given) has the necessary a — > —a symmetry 
and thus contributes to h n . It follows from ( |2.19| ) that the mass of two -Bi-breathers, 
2m 1; is smaller than the one of the soliton-antisoliton pair, 2M, in the region (3 < 1. 
Consequently, 

h n (r) oc exp(-2mir), < (3 < 1 (2.31) 

at large distance r, in agreement with the (3 — > limit ( |2.5| ). It stands to reason that 
subsequently the soliton-antisoliton pair determines the large- distance asymptotic of h n , 



h n {r) oc exp(— 2Mr), 



1 < (3 < 2 



(2.32) 
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where the (3- dependence of M can be deduced from eq. Q2.22|) . At (3 = 2, the result 



( p. 11 ) is recovered. We omit a tedious calculation of the inverse-power law dependence of 



prefactors on distance in (|2.31| ) and (p.32[) 



One concludes that, for a given (3 < 2, the large- distance exponential decay of h n is 
faster than the one of h p . Although the correlation lengths depend continuously on j3 
for both hp and h n , the derivative of the density correlation length with respect to (3 is 
discontinuous at (3 = 1. The correlation lengths coincide just at the collapse (3 = 2 point, 
where, as is clear from the exact formulae ( |2.10| ) and ( [2.11 ), h p and h n differ from one 
another only by inverse-power law prefactors. 

3 Density correlations near a rectilinear hard wall 
3.1 Weak coupling 

Like in Section 2.1, at the lowest order in f3, the direct correlation functions c p and c n are 
given by relations (|2.2|) and fl2.3|) , respectively. However, now, the wavy line K, the sum 
of the chain diagrams, is a function K(x,x'; \y — y'\) which must be calculated with the 
wall taken into account. An arbitrary ew value can be assumed for the dielectric constant 
of the wall material. At the lowest order in (3, K can be computed with the density n(x) 
replaced by its bulk value n. This K has been implicitly considered in refs. [20, 45], where 
it was used as the charge correlation function, h p = K. Except in the cases = oo (ideal 
conductor wall) and ey/ = (ideal dielectric wall), K has an algebraic decay, for large 
\y — y'\, of the form — ew/[^n(y — y') 2 } exp[— n(x + x')]. Therefore, 

U p (x, x'\ \y - y'\) = n 2 K(x, x'\ \y - y'\) ~ €wU exp [-k(x + x')} (3.1) 

7r(y - V ) 

The Meeron diagram (|2.3|) gives for the asymptotic form of the density Ursell function 

U n (x, x'] \y-y'\) = n 2 c n (x, x'\ \y - y'\) ~ - exp [-2k(x + x')] (3.2) 

27r 2 {y -y'Y 
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These formulae can be generalized to v — 2, 3 dimensions: 
U p (x,x';\y\) ~ 



[v 



ewU r / i /\i 

— — — exp [- K (x + x)) 



U n (x,x'; |y|) 
with k = [2(z/- l)n(3n} 1 / 2 . 



-w 



2{v - l)% 2 |y| 



— exp [— 2k(x + x')] 



(3.3) 
(3.4) 



3.2 /?=2 

When (3=2 and ew = 1, the charge and density Ursell functions can be expressed in terms 
of auxiliary functions g a+ , where a = ±, as H 



U p (x,x';\y -y'\) = -2m 2 
U n (x,x';\y -y'\) = 2m 2 



9-+(x,x';y-y')\ + \g++(x, x'\ y - y')\ (3.5a) 



g-+(x, x'; y-y')\ - \g++(x, x'\ y-y 



/\ 1 2 



(3.5b) 



where m = 2irz (z is the fugacity). It will now be shown that, near the wall, each g 
function has a slow algebraic decay for large \y — y'\. 
Each g is a sum of two terms 



g a+ (x, x'- y-y')= g^\v - r') + gZf (x, x<- y - y' 



(3.6) 



The first term in the rhs of (|3.6| ) is the same as in the bulk. It has a fast (exponential) 
decay and does not contribute to the asymptotic form, which is entirely due to the second 
term. This second term is determined by its Fourier transform as 



g w 4\x,x- 



/oo A] 
-g:f(x : x'-l)e W [il(y-y')] 
-oo Z7T 

These Fourier transforms are 



(3.7) 



gjf(x,x';l) 
gjf(x,x';l) 

gZf(x,x';l) 



g-t\x,x';l) 



m 



— exp[— k(x + x')], I < 



m(k — I) 



exp[— k(x + x')}, I > 



2k(k + l) 
k + I 

— — — exp[— k(x + x')], I < 
exp[-k(x + x% I > 



(3.8a) 
(3.8b) 
(3.8c) 
(3.8d) 



16 



where k = (m 2 + Z 2 ) 1 / 2 . Thus, g™+ l (x, x';l) is a function of I singular at I — 0. This 
singularity generates in g™f l (x, x'; y — y') an algebraic decay at large \y — y'\. The corre- 



sponding asymptotic expansion is obtained by splitting the integral in (|3.7| ) into the I < 
and I > contributions, and evaluating each contribution by successive integration per 
partes. For instance, calling for simplicity N(l) the function g™+ l (x, x'\ I) when I < 0, one 
obtains the asymptotic expansion 
r° dl 

/ N (l)exp[il(y-y')} 

= tt" \ —N(0) + t — 1 —^N'(0) + '-—N'UO) + ...) (3.9) 

2vr \ y - y> 1 ; (y - y'Y K ' (y - y'Y K } J 1 ; 

A similar expansion holds for P(l), the function g™+ u (:r, x'\ I) when / > 0. Using these 



expansions finally gives 

\g ++ (x,x';y-y>)f = ^ 2 {j^y 2 + 



1 2(x + x') 
m 2 m 



■1 



x exp[-2m(x + x')} (3.10a) 

and 

\g-+(x,x';y - y')| 2 



2(x + x' 



1 

7 j I + ...jexp[-2m(a; + a; / )] (3.10b) 



47r 2 [ (y — y') 2 m (l/ — u' 
The dominant term in eqs. ( |3.10| ), of order l/(y — y') 2 , determines the asymptotic 
form of the charge Ursell function ( |3.5a|) p0| , 



U p (x,x';\y-y'\)~- - — exp [-2m(x + x')\ (3.11) 

V 7T / (y - y') 2 

But this dominant term cancels out in the expression ( |3.5b| ) of the density Ursell function, 



the asymptotic form of which is governed by the subdominant term in eqs. ( p.l(J| ), of order 
1/ '{y — y'Y ■ The final result is 

U n (x,x'; \y-y'\) ~ 1 —exp[-2m(x + x')] (3.12) 

2tt% - y'Y 

3.3 Conjectures 

The functions U Ptn (r) and S p ^ n (r) differ from one another only by a term containing 
S(r), which has no effect on their identical large- distance behavior. With regard to the 
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definitions of the asymptotic characteristics f p (|1.17|) and f n (|1.25|) , the results of the two 



previous subsections can be summarized, for v = 2 dimensions and ew = 1, by 



n 



f p (x,x') = exp [— k{x + x')\ , (3 — > (3.13a) 



7T 



m N 2 



/ p (x,x') = -( — ) exp [-2m(z + a/)] , (3 = 2 (3.13b) 



and 



/ n (x,x') = ^ exp [-2k(x + x')] , 0->O (3.14a) 

/n(x,x') = exp [-2m(x + a/)] , /? = 2 (3.14b) 

The explicit forms of f p and /„ in the two exactly solvable cases has an appealing 
feature: both charge and density functions factorize in x and x' particle coordinates. This 
can be intuitively explained by the fact that one studies the leading asymptotic y — > oo 
limit of pair correlations along the wall, in which probably there is no correlation between 
x and x' coordinates of particles. In both (3 — > and (3 = 2 cases, the decay of factorized 
functions into the bulk has the exponential form obtained for the corresponding bulk 
functions in subsections 2.1. and 2.2. With regard to the general form-factor analysis in 
subsection 2.3., it is therefore tempting to write down 

2 

f p (x,x') = -—^exp[- mi (x + x')}, 0<(3<2 (3.15) 
with the mass m\ of the elementary I^-breather given by formula (|2.27 ). This formula 



reproduces correctly the two solvable cases ( |3.13j ) and reflects the dominance of the lightest 



particle in the sine-Gordon spectrum in the large-distance behavior of charge correlation 
functions. The value of the prefactor in ([3.15|) is fixed by the sum rule ( |1.18|) with 
v = 2 and ew = 1. As concerns the density asymptotic charasteristic f n (x,x') ( p. 14 ). 



the prefactor of the exponential acquires the same value in the (3 — > limit as well as 
at (3 = 2, what motivates us to suggest, taking into account large-distance asymptotics 
O) and (|232p, that 

f n (x,x') = — - exp [— 2m\{x + x')] , < (3 < 1 (3.16a) 
fn(x,x') = — r exp \-2M{x + x')} , K (3 < 2 (3.16b) 



In particular, at the boundary x = x' = 0, f n (0, 0) is supposed to have a universal value 
l/(27r 2 ) independent of (3. We were not able to give some general argument for such a 
result. 

We emphasize that formulae (|3.15|) and (|3.16|) are only conjectures which must be 
verified, for example, within a systematic weak-coupling expansion in the presence of a 
dielectric wall. Such calculations for the correlation functions are tedious and far from 
simple, and will not be presented here. 



4 Scaling analysis in a disk 

Under the neutrality constraint N + = AL, the grand partition function of the two- 
dimensional symmetric TCP in a domain D, bounded by an impermeable hard wall (for 
simplicity, uncharged and with no image forces, ew = 1) is written as 

oo 2N 

r N 

Qn,n = / n d W 2 n 4 W^\p,n) (4.1b) 



i=i 

where Qn,n is the configuration integral and 

<)(pnH n^| P .-P J |>. ; n J P (42) 

lli,j=l I Pi lv j\ 

denotes the interaction Boltzmann weight of N particles of charge +1 with coordinates 
{Pij^Li and N particles of charge —1 with coordinates {xii\f =1 . From the explicit repre- 
sentations 

oo J2N n N 



1 OO 'IN „ L\ 

N 



x^[5(r- Pl ) + 5(r- ni )] (4.3a) 
i=i 

oo „ N 



1 OO 'IN . JV 

(n(r)n(r')) = § £ ^ \{ d 2 Pi dV <°(P, n) 



N=0 ^ ' V i=l 

A 7 " 

x £ [5(r - Pi ) + 5(r - n,)] £ [5(r' - Pi ) + 5(r' - n*)] (4.3b) 

i=l j=l 
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etc., one readily gets the important relations 

d 

d 2 r n(r) = z— InH (4.4a) 

d oz 

[ dVS„(r,r') = z— n(r) (4.4b) 

etc. 

The domain of interest in this section is the disk of radius R, D = {|r| < R}. When 
one rescales the particle coordinates as = Rp' { and rij = Rn 1 ^ it becomes evident that 
R) depends only on the dimensionless combination z 2 R 4 ~ 13 . As a consequence, 

*§-('-!)■§ 

If one accepts the expected value of c = — 1 in the universal lni?-term of the large- R 
expansion ( 1.24 ) (% = 1 for the disk), InH = — /3fi takes the form 



In E(z, [3, R) = (nR 2 )pp - (2-kR)^ - - In f z^^-Wr) + const + • • • (4.6) 

where 

PP = fv(P)z 1/(1 - m (4.7) 

p is the bulk pressure, and 

/? 7 = f s {P)z 1/{2 - m (4.8) 
7 is the surface tension. The bulk particle density n is given by the equation 

n = z dm = i (49) 

which gives the equation of state ( |1.11|) . 

The particle density n(r, R) depends on both the radius R and the distance < r < R 
from the centre of the disk. By using the above mentioned scaling transformation of 
particle coordinates, one has 



n 



(r, R) = zW-V^g (z x 'W>r, z^'^r) (4.10) 
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with an unknown function g. The origin can be moved to the boundary via the coordinate 
transformation x = R — r. In order to distinguish between the two different functions, we 
will use the obvious notation 

n R {x) = n(R-x,R) (4.11) 

Similarly as in ( |4.10| ), 



n R (x) = zW-fMg (z^-^x, z^^R) (4.12) 

with an unknown function g different from g. The definition ( [4.1 1|) implies 

dn(r,R) dn R (x) 
dr dx 



(4.13a) 



dnjr, R) dn R (x) dn R (x) 

8R 8x 8R 1 } 

The transition from the disk of radius R to a rectilinear hard wall can be understood as 
the limiting R — > oo procedure, with liniR_ +0O n R (x) = n(x) being the particle density at 
a finite distance x > from the plain hard wall. 

It is evident that the derivative of nili d 2 pjd 2 nj Wp (p, n) with respect to R is 
2vri? /<f nili d 2 p<d 2 7i i W / S JV) (p, n) E^i [<J(Pi - R) + <J(n< - R)] where R is a position vec- 
tor of any point at the disk boundary. Hereinafter, the integration range from to R 
formally means that < {|pj|, |nj|} < R. As a result, 

^^5 = 2^71^(0) (4.14) 



With regard to eqs. (|4.6j ) and (|4.9|) , the finite-size correction of the particle density at the 
boundary reads 

»«("> = f 1 -!)»- 1<*> - lab + - (415) 



The well-known contact theorem n(0) = [1 — (j3/A)\n [46, 47 1 results as the R — > oo limit 
ofeq. flPgp . 



By using the same procedure for H n(r, i?) one arrives at 

^ r ~ n(r, i?)] = ER ^ dip [(n(r)h(R, <p)) - n(r)8(r - R)] (4.16) 

J — n 



dR 
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where, in polar coordinates, R = (R, <p). Using eq. (|4.13b| ) and after some simple algebra, 
one can pass from Q4.16Q to 

d -^ + ^-«J>^ (4,7) 

In the limit R — > oo and for a fixed x, ur{x) = n(x)+0(l/R), so lim^oo dnn(x) / dR — > 0. 
When one introduces the y-coordinate as follows y = Rip, eq. ( j4.171) takes the form 

^ = r dyU n {0,x;y) (4.18) 

OX J-oo 

and we recover the two-dimensional version of the WLMB equation [IS 



We are now ready to take advantage of the exact grand-canonical relation (|4.4b ). With 
regard to the scaling form of ur(x), eq. (|4.12j ), it holds 



In the limit R —>■ oo, the last term on the rhs of eq. (|4.19|) vanishes. Inserting the resulting 



zdn(x)/dz into ( 4.4b ), one finally obtains 



OO POO 



i / / i n i i \ n ( x ) x dnix) .. 

. da// dy S n {x,x'\y) = ttttt + —, ttt, » » ( 4 - 20 ) 

o J-oo V ' l-(/9/4) 2[1 - (/?/4)] dx v ; 

This is the local form of the compressibility sum rule in the presence of a plain hard wall. 
Indeed, in the limit x — > oo, n(x) — ■> n and dn(x)/dx — > faster than any power-law due 
to screening. Consequently, in the bulk, (|4.20D reproduces (|1.22| ). It is straightforward to 
prove that eq. (|4.20|) is valid for any value of eyy- 

5 Mobius invariance and sum rules 

We now consider the TCP confined to a two-dimensional domain D, and study the action 
of the Mobius conformal transformation 

, az + b dz' — b . . 

z = — —j) z= 7— ( 5A ) 

cz + a —cz' + a 
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(with free complex parameters ad — be 7^ 0) of complex particle coordinates (z, z) on the 
configuration integral Q n,n Q4-irj| ) . Under the conformal transformation ( |5.1|) , the domain 
D is mapped onto the one denoted as D', the surface element dzdz is written as 

Azi _ = ^d-bc)(U-lc) 
(a — cz') z (a — cz'Y 



the square of the distance between two particles takes the form 

u z .i2 = (ad-bc)(ad-bc) 2 

1 3 [a — cz'j){a — cz'i){a — cz'j)(a — cz'a) 1 3 

so that the Boltzmann factor ( (4.2[) reads 

-N/3/2 



W^ N \p,n) = [(ad-bc)(ad-bc) 



N 

x II [(« - aMfi ~ rt)f 2 [(a - Oifl ~ Of 2 <V, n') (5.4) 



i=i 

Finally, putting c — c—1 and a = r , a = f , one finds 
r N 

J D I[d 2 p t d 2 n t W { p N \p,n) = [(b-r d)(b-f d) 



N[2-((3/2)] 



i=l 

X 



" d 2 Pt d 2 n t {N) 



This means that the configuration integral of the TCP Coulomb system confined to the 
domain D is mappable onto the one of the TCP Coulomb system (with the same number 
of charged particles and at the same inverse temperature (3) confined to the new domain 
D' and being under an additional action of an external gravitational source with specific 
coupling strength equal to 4 — (3. Since in the stability range of (3 < 2 the gravitational 
coupling 4 — (3 > 2 and the lhs of eq. ( |5.5| ) is finite, the formalism must ensure that, 
for any values of free complex parameters b and d such that (b/d) ^ r , the gravitational 
source is localized outside of the domain D' in order to prevent the divergence of the rhs of 
eq. (|5.5| ) due to the gravitational collapse. According to the definition ( |4.1a|) , the relation 



between the configuration integrals ( [5.5|) also implies an analogous relation between the 
grand partition functions, 

E(z,(3,D) = E(z,(3,D'\r ) (5.6a) 
(b-r d)(b-f d)] lH/3/4) (5.6b) 
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Taking for the domain a disk with its center at the origin, D — {|r| < R}, let us set 
b = R 2 and d = f (r f 7^ -R 2 ), besides the already chosen c = 1 and a = r , so that the 
Mobius conformal transformation (15.11) takes the form 



, r z + R 2 r z' - R 2 

z = — , z = (5.7) 

z + r —z' + r 

It is straightforward to show that under this transformation the disk zz < R 2 maps onto 
a "dual" domain D' defined by the inequality 

(R 2 -r f )(R 2 -zz') < (5.8) 

If the gravitational point is outside of the disk, i.e. r f > R 2 , the transformation (|5.7|) 
maps the disk onto itself. Relation (|5.5|) results into 



rR " 

jf n4d 2 « i < ) (p,n) = (r f -R 2 ) 



i=i 

x 



R TT ^ ^ W (N) (x> n) (5 9) 

If tqTq < R 2 , the disk is mapped onto its complement in the two-dimensional space and 



eq. (|5.5|) gives 

N 



d 2 pi d 2 rii rrAN). 



J ndWn^ffp,!!) = (i^-/-„/-n 



i=i 

00 A ,1-,, ,1-' 

X ' 



f°° TT Cl ft CIJI, T zr(JV)/ V r in"\ 



'« i=i l r o - I |r - n, 
In both cases, the gravitational source lies outside of the dual domain, as is required 
by the condition of stability (see the comment in the above paragraph). There exists a 
special Mobius transformation which maps the disk directly onto a half-plane, but we will 
not discuss this case. 

Although the present paragraph will not be used in the following, let us remark that 
the self-duality of the disk system ( |5.9| ) produces a relation between the grand partition 
functions of type (|5.6|), with D' = D = disk and z = z(r f — R 2 ) 2 ~^^ 2 \ Since S(z, ft, R) 
depends on ft and the combination z 2 R 4 ~ 13 , E(z, ft, R\r ) will depend on ft and the combi- 
nation z 2 R^~P I {tqTq — R 2 ) A ~P . Let x = \ro\ — R be the distance of the gravitational point 
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from the surface of the disk. The consequent ratio R/(2xR + x 2 ) diverges in the limit 
R — > oo, x — > 0. In this case, the asymptotic formula (|4.6|) can be applied. When x is 
finite, the asymptotic formula fl4.6|) is applicable only when one multiplies z by an appro- 
priate independent constant via changing the zero reference energy of the gravitational 
potential. 

We will now use the self-dual relation ( [5.9D to derive a new sum rule Q5.14| ) and its 
more detailed version ( |5.17| ) for the density structure function S n in the case of the disk 
geometry. First, we rewrite the rhs of ( |5.9| ) as follows 



r n 



R 2 /(r f ) 



Pi/r )(l -Pi/r ) 



2-09/2) 



d 2 rii 



R 2 /(r fo) 



ni/r )(l - rii/ro) 



2-08/2) 



(5.11) 



Then, since R 2 /(r r ) < 1, \pi/r \ < 1 and \rii/r \ < 1, we perform the large-|r | expansion 



1 - R 2 /(r f ) 
- z/r )(l - z/f ) 



2-09/2) 



-UA— 

\ 2 / r r 



2-f 



zz 



r r 



+ 



(5.12) 



for each z = {pi,rii}. The term of order (rofo) in ( ^.11| ) exactly reproduces the lhs of 
( |5.9p , and the coefficients of higher-order terms in inverse powers of r and f must be 
identically equal to zero. The terms l/r and l/f trivially vanish. The l/(r f ) term 
vanishes when 



rR n N 

/ J] d 2 Pi d 2 ni (Pj + n i)(Pk + n k )W^ N) (p, n) 

J ° 8=1 j,k = l 



R A 



^/;ndW<>(P,) 

(5.13) 



Within the grand-canonical formalism ( (4.1|) -( ^3|) , the equality ( p,13| ) is equivalent to 



/ d 2 r / dVrT^(r,r') 
Jo jo 



<9 

-z— lnE(z,P,R) 



(5.14) 



2-(P/2) dz 

Everything that has been done in the case of the configuration integral can be adapted 

to the particle density. Let us "rotate" the Mobius transformation ( |5.7|) around the origin 

by the multiplication factor (fo/ro) in order to obtain the identity z' = z in the limit 

|r | — > oo, and denote by a the ratio R/r whose absolute value is smaller than 1, 

. z + aR 



1 + (az/R) 



(5.15) 
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For a given point r = (r, r) in the interior of the disk, the previously developed formalism 
implies 



f n d W £ [5(v - Pj ) + 5(v - n,)] Wf\p, n) 

J0 1=1 7=1 



'1 - aa) 2 



1 + ar/R) 2 (l + af/R) 2 Jo 
1 — aa 



r n 
IK* 

i=l 
2-09/2) JV 



1 — aa 



1 - api/R)(l - api/i?) 



2-08/2) 



d n, : 



E " Pi) + 5(r - %)] W^V n) (5.16) 



^1 — arii/ R){1 — arii/R) 

When one expands the rhs of Q5.16D in powers of a = a x + ia y , \a\ < 1, the requirement of 
the nullity of the coefficients attached to a x and a y leads to the relation 

(4 - (3) [ R d 2 r' r • r'S^r, r') = 4r 2 n(r, R) + (r 2 - R 2 )r^-n(r, R) (5.17) 
Jo ar 

The new sum rule ( |5.14| ) results from the more detailed eq. ( |5.17[ ) when one integrates 
J R d 2 r both sides of the latter, then uses an integration per partes and finally applies the 
equality (|4.4a| ) . 

This new sum rule (|5.17|) has a consequence for the density near the boundary. Let 
us divide eq. ( |5.17| ) by |r|, and move the origin to the boundary via the coordinate 
transformation r = R — x described in section 4. Then, using eq. (|4.19| ), and in the 
R — > oo limit, we find after some simple algebra that 



n R(x) — n ( x ) ^ { (4 — f3) / dx' x' / dy S n (x, x'\ y) — Axn(x) — x 

2R Jo J-oo 



,dn(x) 
dx 



+ ■■■ 

(5.18) 



This formula determines the leading correction of the density tir(x), at a fixed distance 
x from the wall, with respect to its asymptotic R — > oo value n(x) due to the curvature 
of the confining disk domain. The quantities on the rhs of ( |5.18| ) are the ones evaluated 
for the rectilinear hard wall. 

The asymptotic formula (|5.18| ) can be explicitly checked at the boundary, where it 
gives 

roo roo 

/ dxx dyS n (0,x; y) + ■ ■ ■ (5.19) 

JO J-oo 



n R (0) = n(0) 



A-(3 



2R 
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Multiplying eq. ( 4.18|) by x, integrating then over x from to oo and performing an 
integration per partes, one finds 



dx [n(x) — n] = dx x dy [S n (0, x; y) — n(x)8(x)5(y)] 



(5.20) 



On the other hand, according to the relation ( 4.4a|) , it holds for the disk 



1 r / \ n 111 t— i q 

d r n(r, R) — n\ = z — nirR 

o dz 



(5.21) 



Assuming the large-i? expansion (H) and passing to the boundary- distance variable 
x = R — r, eq. ( |5.21| ) takes the form 

r« 9(/3 7 ) 1 



2tt / dx(i? - x) \n R (x) - n] = -(2ttR)z- 
Jo 



+ ■ 



dz 12(1 - (/3/4)) 
In the limit R oo, this equation implies the obvious boundary relation 



(5.22) 



/•oo ( 

/ dx [n(x) — n] — —z- 



dz 



(5.23) 



Taking into account the z-dependence of /?7 ( |4.8| ) in ( p.23| ) and considering the previously 
derived relations ( |5.19| ) and ( p.20| ), large- R asymptotics ( |4.15| ) is reproduced correctly up 
to the (l/_R)-term, which confirms the validity of ( |5.19| ). 

The asymptotic formula ( |5.18| ) implies a generalization of the second moment den- 



sity sum rule ( 1.23 ) to the case of a rectilinear wall. Indeed, the local version of the 
compressibility sum rule ( [4.20D can be used for rewriting Q5.18D as 



n R (x) = n(x)-^^(A-(3)J o C °dx'(x'-x) J°^dy S n (x,x f ;y) + x 2 ^^-^ + ■ ■ ■ (5.24) 

We will suppose that the half-infinite Coulomb system has good screening properties into 
the bulk, i.e., [n(x) — n] decays faster than any inverse-power law as x — > oo and all 
moments dx x l [n(x) —n] exist. From eqs. ( |5.22| ) and ( |5.23| ), one then gets in the limit 
of large R 



/•OO /"OO 

2ttR / dx [nji(x) — n{x)} — 2n dx x [nix) — n] 
Jo Jo 



12(1 - (0/4)) 



(5.25) 
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Inserting ( |5.24|) into (|5.25|) and after some algebra we arrive at 

roo roo roo j 

Jo dx L dAx ' -x)]_JySn{^y) = 48?r(1 _ m))2 (5-26) 

Actually, it can be shown that the boundary sum rule ( |5.26| ) can be derived directly 
from the second-moment sum rule in the bulk ( |1.23| ) by using a simple assumption. Indeed, 
S n (x, x'\ y) can be decomposed as the sum of the bulk structure factor plus a surface term: 

S n (x, x'; y) = 5 n bulk (r) + S^ix, x'; y) (5.27) 

where r = [(x — x') 2 + y 2 ] 1 ^ 2 and S^ ulk (r) is the bulk density structure factor appearing in 
( |1.23| ) [ (|5.27 ) is just a definition of e ]. When this decomposition is used in the lhs of 



( |5.26| ), assuming that S 1 ™ e (x, x'; y) has a fast decay when x or x' or both go to infinity 
(see for instance the explicit expressions in the weak-coupling or (3 = 2 cases considered 
in sections 3.1 and 3.2), the corresponding integral is absolutely convergent, the order of 
the integrations on x and x' can be freely interchanged, and since (x — x')S^ irface (x, x'\ y) 
is odd under the interchange of x and x', this integral vanishes. This reasoning however 
does not apply to the contribution of the bulk structure factor which does not decay to 
when both x and x' go to infinity for a fixed value of x — x'; the corresponding integral 
is not absolutely convergent, and the order of the integrations cannot be changed. Thus, 
only the bulk part of S n contributes to the lhs / of (|5.26|) which can be rewritten as 



POO POO POO 

1= dx dss dy S™ k (\s\;y) (5.28) 

^0 J —x J —oo 

where we have used the integration variable s = x' — x rather than x' and made explicit 
that S^ ulk depends on x and x' only through \s\. Performing the integration on x per 
partes gives 



/oo roo x— oo ^oo poo 

dss dy S^ lk (\s\-y) + dx dyx 2 S™ k (\x\; y) (5.29) 
-x J — co x=0 JO J —oo 



Because S% has a fast (exponential) decay at infinity, the first term in the rhs of ( |5.29| ) 
vanishes. Since S^ ulk (\x\;y) is a function of r = (x 2 + y 2 ) 1 ^ 2 only, / is 1/4 of the lhs of 
( |1.23j ), and (|1.23|) results into ( |5.26|) . It should be noted that the above derivation of the 
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boundary sum rule (|5.26| ) still holds for an arbitrary value of the wall dielectric constant 
ew, including the special cases ew = and ew = oo. 



6 Universal finite-size correction for the disk 



In most considerations of the previous two sections, we have assumed the validity of the 
large- R expansion of InH ( (4.6p for the disk, in particular the value —1/6 of the coefficient 
of the universal In i?-term. Here, we prove this assumption. 
Using the formula |r — r'| 2 = |r| 2 + |r'| 2 — 2r • r', we write 

"R - R 2 



R d 2 r f R d 2 r'\r-r'\ 2 Si R \r,r' 
o Jo 



Z dz 



I d 2 r\r\ 2 n(r, R) — — In 5 

Jo 1 1 V ; 2 - ((3/2) 



(6.1) 



where we have used the relation ( |4.4b| ) and the new sum rule ( |5.14| ). Let us consider the 
large- R expansion (|4.6|) in a general form 



InS = (nR 2 ) (l-^\ n - (27ri2)/3 7 + / 



(6.2) 



where n oc z 1 ^ 1 /3j oc z 1 ^ 2 ^l 2 ^ and / is an as-yet-undetermined function of (3 and 
the combination z 2 R 4 ~ 13 . The representation ( |6.2|) when combined with (|5.21|) implies 

1 - (P/A) I df 



Pi 



nR 



i z ^~ ~ J d 2 r [n(r, R) — n] j 



Substituting (3^( in ( |6.2| ) by the expression ( |6.3| ) and then inserting InS into (|Q 
obtains 



(6.3) 



one 



Td 2 r TdVlr-rf^^r') 
tcR 2 Jo Jo 



2 d f R ( r 2 \ 
-n Z dzJo dV Hr,R)-n) 



d . 2 f d 
-z—f + -\z- 



tt(1 - ((3 /A)ydz J ' vr \~dz) f (6 ' 4) 
This equation, which is exact for an arbitrary disk radius R, will now be analyzed in the 
limit R — > oo. In that limit, since S^'(r,r') decays along the boundary like l/|r — r'| 4 
(see section 3), there is no relevant surface contribution to the integral on the lhs of ( |6.4|) . 
and we can make use of the bulk sum rule ( |1.23|) derived in [fZ5] , 



lim At f R d 2 r [* dVlr - r'| 2 ^(r, r') = — - 

fi-oo nR 2 Jo Jo 1 12vr 1 - 



(1 - (/3/4)) 2 



(6.5) 
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As concerns the rhs of ( |6.4|) , the unknown integral over the disk area can be represented 
in the x = R — r coordinate as 

J Q R dV ( X " ^) K r > R)~n] = 4vr J* dx x (l - |) (l + ^) [n R (x) - n] (6.6) 

In the limit R —>■ oo, due to good screening properties of the plasma, only the term 
oc / °° dx x[n(x) — n] survives. Based on the scaling analysis of section 3, since n(x) = 
n 9{Vn x )i this integral is dimensionless, and therefore its derivative with respect to z 
vanishes. Consequently, without the integral on the rhs of eq. ( |6.4j ), the solution of that 
equation combined with (|6.5| ) reads 

f = --\n (Rz^-W) + const + • • • (6.7) 

in full agreement with the expected universal large- R behavior ( |4.6| ). 



7 Concluding remarks 

The pecularities of the charge correlation in a Coulomb fluid are determined by the screen- 
ing effect caused by the long-ranged tail of the Coulomb potential. The density correlation 
function in a charge-symmetrical two-component plasma, decoupled from the charge one 
at the level of the OZ equation [see relations ([2.11) 1, does not feel this tail directly. Yet, 
the density correlation bears some similarity with the charge correlation. Although the 
density correlation decays always faster than the charge one, it still exhibits an alge- 
braically slow decay ~ |y| _2iy ( |3.4|) at large distances |y| along a rectilinear hard wall in 
any dimension v. An analogous phenomenon for the charge correlation, which decays like 



~ \y\~ v ( |3.3|) , is related to an asymmetry of the screening cloud near the wall. 

In this paper, we have concentrated mainly on the two-dimensional TCP of pointlike 
charges. The form-factor analysis of the equivalent sine-Gordon model was done in the 
bulk case. In contrast to the charge correlation, the density one exhibits a discontinuity 
of the slope of the correlation length, namely at point (3 = 1 [see relations ( |2.31|) and 



( p. 32] )]. At the collapse point (3 = 2, as seen in eqs. Q2.10D and fl2.11|) , the charge and 
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density correlation lengths become the same and the difference in the asymptotic decay 
is only in the inverse-power law prefactors. 

The TCP in contact with a rectilinear hard wall of dielectric constant ew was mapped 
onto an integrable boundary sine-Gordon model in two cases: ew — > oo (ideal metal wall, 
Dirichlet boundary condition) |J and ew = (ideal dielectric wall, Neumann boundary 
condition) |5]]. Interestingly, just in these two cases the charge and density correlation 
functions decay exponentially also along the boundary. The next potentially integrable 
case corresponds to the plain hard wall with ew = 1: at least, the explicit solution is 
available at the free-fermion point 0. For this case, we have suggested an explicit form 
of the asymptotic characteristics along the wall, f p (x,x') (|3.15|) and f n (x,x') ( |3.16| ), as 
a natural interpolation between the exact results in the f3 — > limit and at /3 — 2. The 
universal value of the prefactor = l/(2n 2 ) in f n deserves attention. To check the validity 
of our conjectures, one should go beyond the Debye-Hiickel limit. 

The scaling analysis of section 4 leads to a new local version of the compressibility 
sum rule in the presence of a hard wall ( f4.2(J| ). 

In the crucial section 5, a symmetry of the TCP confined to a disk with respect to the 
Mobius conformal transformation of particle coordinates, which induces a gravitational 
source point interacting with the particles, is established. This symmetry is of interest by 
itself: it produces a path between the pure Coulomb system and the one under the action 
of a gravitational point. The gravitational point acts in the same way on positive and 
negative charges, so that the total particle density is relevant. As a result, the new sum 
rule (|5.14 ) for the density correlation function and its local version Q5.17 ) are derived for 
the disk geometry of the Coulomb system. The finite-size formula (|5.18| ) for the particle 
density at a fixed distance from the wall and the boundary version fl5.26| ) of the bulk 
second moment of S n ( |1.23[ ) also deserve attention. 

The new density sum rule (|5.14|) is used in section 6 to derive the universal finite-size 
correction term of the grand potential for the disk geometry. We would like to stress 
a fundamental difference between the derivation of the universal finite-size correction 
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term for a system which is finite in each direction (a disk or a sphere, a density sum 
rule is relevant) and for a system which is infinite at least in one direction (a strip or a 
v- dimensional slab, a charge sum rule is relevant |33| Q). 

In most cases, we have restricted ourselves to the symmetrical-charge TCP and ew = 1 
in order to maintain the clarity of presentation. The extension of the formalism to the 
asymmetric TCP and arbitrary ey/ is usually straightforward. 
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